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We show how to obtain the ultraviolet central charge from the exact S'-matrix for a 
wide variety of models with a U(l) symmetry. This is done by coupling the U(l) current 
J to a background field. In an N=2 superconformal theory with J the fermion number 
current, the OPE of J with itself and hence the free energy are proportional to c. By 
deforming the super symmetry into affine SM(2) g quantum-group symmetry, this result can 
be generalized to many £7(l)-invariant theories, including the iV=0 and N—l sine-Gordon 
models and the SU (2)k WZW models. This provides a consistency check on a conjectured 
S'-matrix completely independent of the finite-size effects expressed in terms of dilogarithms 
resulting from the thermodynamic Bethe ansatz. 

7/93 



' address after Sep. 1, 1993: Department of Physics, University of Southern California, Los 
Angeles, CA 90089 



One way of solving integrable two-dimensional field theories is to put the model on 
the lattice and diagonalize the quantum Hamiltonian by using the Bethe ansatz, much 
in the same way that one solves a quantum spin chain. However, another approach has 
become popular in recent years: one uses the enormous constraints of an integrable theory 
to conjecture the exact S-matrix [Q. Instead of dealing with field theories, one studies 
relativistic quantum particles. From the exact S-matrix, one can derive exact form factors 
and extract correlators. One can also derive the free energy at non-zero temperature 
from the thermodynamic Bethe ansatz (TBA) ||. An advantage of the exact S-matrix 
approach (and the main reason so much effort has been devoted to it in recent years) is 
that models can be solved by defining them as conformal field theories perturbed away 
from the critical point, without knowing the action or Hamiltonian explicitly 

When conjecturing an exact S-matrix, it is important to check that it has the 
right properties by comparing results with conventional Lagrangian perturbation theory, 
conformal perturbation theory or numerical simulations. One very unambiguous check 
is to calculate the Casimir energy from the TBA (by changing the roles of space and 
time, finite-temperature effects become finite-size effects), because in a critical theory, it 
is proportional to the central charge c || . Many such calculations have been done by now 
(too many to refer to them all). The result is given in terms of sums of dilogarithms, and 
to obtain the final result, one must use or derive some marvelous dilogarithm identities. 

The purpose of this letter is to show that there is another simple quantity which 
is proportional to the central charge in any theory with N=2 supersymmetry, and many 
theories with a U(l) symmetry. We show how to calculate this quantity from the S'-matrix, 
giving a check on the S'-matrix. This calculation is done at zero temperature but requires 
adding a background field coupled to a conserved charge. 

In the past, the effects of a background field coupled to a U(l) current have been 
used to relate parameters in the Lagrangian to the S- matrix H0J§. Here we will see how 
there is an even bigger payoff in theories with N=2 supersymmetry, where the U(l) is 
the fermion number. Because the chiral anomaly in an N—2 theory is proportional to 
the central charge c, c is determined directly from this calculation, independently of any 
details of the specific Lagrangian. For example, no renormalization effects need be studied. 
Moreover, this is not the only case where this observation is useful, because one can deform 
the N=2 supersymmetry into an affine su(2) q quantum-group symmetry while keeping the 
model integrable 0. (If one bosonizes the U(l) current, this is equivalent to changing the 
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radius of the boson.) Many integrable models with a U(l) symmetry can be obtained in 
this manner, so this provides a strong check on these S'-matrices as well. 

Consider a theory with a U(l) symmetry generated by a current J^. We write the 
conservation law as OJl + OJr = 0. The corresponding conserved charge Q can be coupled 
to a background "electric" field A which we take to be constant to ensure that we don't 
destroy the integrability of our theory. This modifies the Hamiltonian to be H = Hq + QA. 
Because A has dimensions of mass, the strength of A controls the position of our theory on 
its renormalization group trajectory and, in particular, in the limit of large A the theory 
is driven to the UV. The ground-state energy of a theory in the limit of large background 
field A thus is related to a property of its UV conformal fixed point theory. 

We first show that the ground-state energy in a background field is related to the 
chiral anomaly of the U(l) current. Suppose at a conformal fixed point of the flow that 
the left and right charges are separately conserved, i.e. <9Jl=<9Jjj=0. We can then couple 
independent background fields to the left and right charges, modifying the hamiltonian as 
H = Hq + AlQl + ArQr. Off the fixed point we will need to set Al = Ar because only 
Q = Ql + Qr is conserved, but for the moment we can treat them independently. The 
action is modified by the background fields as 

S = S + [ ^(A L J L + ArJr). (1) 

J ^7T 

We regulate the theory by taking spacetime to be a torus of euclidean time (3 and 
spatial length L; at the end we take the volume to infinity. In the large (3 limit we 
have Tre~@ H — *> e -P E ( A ) ; where E(A) is the ground-state energy in background field A. 
Evaluating the partition function as a path integral using (P , we find in this limit 

e -P(E(A)-E(0)) = u- J ^(AlJl+ArJr)^ ^ 

where the correlator is evaluated in the theory without the background fields, i.e. using 
Sq. By Wick's theorem, (^) yields 

E(A)-E(0) = ± J^J ^((A l J l (z)+ArJr(z))(A l J l (w)+A r Jr(w))}. (3) 

The operator products of the currents appearing in (|3|) are generally of the form 

(J L (z)J L (w)) = - — + (Jr(z)J r (w)) = K + 

(z — w) z V {z — w) z V 

(Jl(z)Jr(w)} = 2ttk5 2 (z -w) — , 



where k is a constant measuring the strength of the chiral anomaly and V = f3L is the 
area of the worldsheet. These OPEs require some explanation. In the infinite volume 
limit we can use the planar expressions for the OPEs, which are the first terms in the 
above equations. The contact term in (JlJr) is required for vector current conservation 
dJL + dJn = 0, for example in ((<9Jl + 9Jr) Jl). The constant 1/V terms in reflect the 
need to absorb the zero mode associated with the normalizable constant modes existing 
in finite volume. The above expressions can be obtained using bosonization by taking 
derivatives of the propagator of a free boson on the torus. From the explicit form of this 
propagator (for example, see Jl0|]) we find two terms, one involving the theta function 9\ 
which results in the first terms in the above expressions and another (necessary for modular 
invariance) yielding the 1/V terms. 

In our infinite volume limit one might be tempted to drop the 1/V terms in (§). 
Plugging (f|) into (|3|) we see that we can do so if the background fields go to zero at infinity 
sufficiently rapidly so that V~ x J d 2 zA(z) — > 0. The remaining pieces yield the familiar JTT] 



result k f (d 2 z /2rr)F (<9<9) _1 F , where F = 8Al — 8Ar is the gauge-invariant field strength. 
In this expression left and right movers are coupled through a AlAr term coming from the 
8 2 (z — w) in (f|). This is the familar statement of the chiral anomaly: by imposing vector 
current conservation we ensure gauge invariance at the expense of left-right decoupling. 
This is why k measures the chiral anomaly. 

We, however, are interested in constant background fields so we clearly need to keep 
the 1/V terms in (|||). The l/(z — w) 2 and 1/(1 — w) 2 pieces vanish in the integrals and we 
are left with the final result for the energy density due to the constant background fields 

elA)s m^m = .^ A i +Ai> . (5) 

The fact that the AlAr term dropped out is sensible: the chiral anomaly is proportional 
to the derivatives of the background fields so for truly constant fields it vanishes and the 
left and right movers should decouple. To make contact with our theory away from the 
critical point we set Al=Ar=A and (|D becomes 

S(A^oo) = -^A 2 , (6) 

where k uv is the chiral anomaly of the UV critical theory as defined in (^). For large but 
finite A, (0) will have subleading terms in (M/A), where M is the mass scale. 



In order to make the connection between k and the central charge of the UV fixed point 
we note that in many integrable models the U(l) symmetry is part of a su{2) q symmetry 
(N=2 supersymmetry is a special case of the quantum-group symmetry at q 2 = — 1.). Our 
convention will be to always normalize the U(l) charge so that the generators of su(2) q 
have charge ±1. We can bosonize the U(l) symmetry as: 

where (j) ~ (j) + 2nR and the normalization constant A is independent of R and determined 
by our above convention. From this we find k = (X/2tcR) 2 . In an N=2 theory, the 
supersymmetry algebra requires that k = c/3. Now the crucial point is that we can change 
the theory by changing the value of the boson's radius R; this changes q. Because changing 
the radius R doesn't change the central charge, this yields a general relation between k 
and c, namely 

S{A^oo) = -^^fA\ (8) 

We now show how to calculate the ground-state energy (H) directly from the S'-matrix, 
giving the promised check. Turning on the background field has the effect of changing the 
ground state. If one were working in field theory or in a lattice model, this would change 
the Dirac or Fermi sea. In our exact particle description, a similar thing happens: the 
ground state no longer is the empty state — it has a sea of real particles. We use rapidity 
variables so that a particle of charge q and mass m has energy qA + m cosh 9. When qA is 
negative, for small enough 9 this particle has negative energy and it will lower the energy 
as part of the ground state. In all known integrable models save that of a free boson, the 
particles behave as if they were fermions; i.e. there is only one particle to a level. Thus 
in the ground state such particles fill all levels with \9\ less than some number B (if the 
particles did not interact, we would have \qA\ = m cosh 5). Any type of particle with 
q a A < can appear in the ground state, so we label them by the index a. Defining the 
associated ground-state densities by p a (9), the energy density of the ground state in a box 
of length L is thus 

£(A) = jYl I a d9 p a (9)(-\q a A\+m a coshe). (9) 

In a general interacting theory, one must resort to perturbation theory. However, the 
exact S'-matrix allows us to derive equations for the p a (9) and B a . These follow simply 
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from putting the system in a box of length L and imposing periodic boundary conditions; 
similar calculations have been done in lattice models ]1^,|13],[14| . It is possible to study this 



constraint because the scattering is completely elastic; the individual momenta remain the 
same in a collision. The scattering need not be diagonal in general — internal labels can 
change. When the scattering is diagonal among particles in the ground state (i.e. the only 
scattering is of the type 0(6*1)6(^2) — > S ab (0\ — 02)b(^)a(9i)), the momentum of the a 
particle of type is quantized via 



gimfc sinh 8 k L 



l[s aka .(e k -e j ) = i, (10) 



an interacting-model generalization of the one-particle relation = In the large 

L limit, taking the derivative of the log of ( |TUD gives the distribution P a (9) of allowed 
rapidity levels for particles of species a. Since we are at zero temperature, the particles fill 
all allowed levels up to B a , i.e. p a (9) = P a (9) for 16*1 < B a , giving 

Pa(0) = ^COsh# + V dd' Pb {d')cj> ah {d-d'l (11) 

2tt V J-B n 



where 



i dlnS ab (9) 

4>ab{V) - 



2tt 89 

The equations fllTD and @ and the equations for B a (given by minimizing the energy with 
respect to B a ) can be put in a more convenient form by defining the "dressed" particle 
energies e a {9) as 

e a (6) = \q a A\ -m a cosh 6 + J2 / dd' <f> ab (9 - 9')e b (e') . (12) 

Substituting this into (||) and using (^T[) one finds that 

£(A) = -Y— [ a d9 cosh6» e a {9). (13) 

It is obvious that the B a are determined by the boundary conditions e a (±B a ) = 0. As we 
will discuss below, the same equations are valid even when the scattering among particles in 
the ground state is not diagonal if the appropriate zero-mass, zero-charge "pseudoparticles" 
are included. 
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The relations (]T^)-(|T3|) are not solvable in general, but the results can be written 
as a series in m a /A by using a generalized Weiner-Hopf technique ||14|| . However, in the 
conformal limit where we make contact with the field dependence derived in (^|), they can 
be solved. Following M, we define e a by 



B 



b 



e a (9) = \q a A\ - -V + J2 / M'MO - d')e b (9'). (14) 
1 b J -°° 

To reach the ultraviolet limit, we take A very large. In this case, the B a are very large 
and the energy (|13[) can be written as 

1^ d9 e e e a {9)- (15) 

^ 2lT J-oo 

Taking the derivative of flT4] ) with respect to 9 and using it to get rid of the m a exp 9 in 
(|15D, after using ( |14"D again one finds 

£(A^^) = -J2 l - S ^U-oo) (16) 

a 

It follows from ( |T4] ) that the e a (— oo) follow from the matrix equation 

(Sab - N ab )e b (-oo) = \q a A\. (17) 

where we define A^ a ;, = J^° oQ d9(f) ab (9). We thus have arrived at a very simple answer for 
the ground state energy at the ultraviolet fixed point. It is given very simply in term of 
the charges and the iV a 5, and the latter follow simply from the ^-matrix. The ground-state 



energy (16) must of course must agree with (||), giving a simple check on the S'-matrix. 

We devote the remainder of this letter to providing examples. First we discuss the 
sine-Gordon model, 

S = J d 2 z [\{d4>) 2 + ^ 2 cos/3(/>] , (18) 

repeating the analysis of || . The conserved charge is proportional to the topological charge 
and the corresponding currents may be written as (0) where the radius is that of the SG 
boson R = f3~ 1 . For (3 2 > 4n the entire spectrum consists of a soliton and antisoliton. 
There is a su(2) q quantum group symmetry with q = — exp(— iiz/t) where the parameter 
t, which appears in the S-matrix element written below, is related to (3 by 
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In our normalization where the quantum-group generators have charge ±1, the soliton 
has charge |, while the antisoliton has charge — \ 0. The N=2 value of (3 (t=2) is 
(3 2 = (2/3)87r; there the charge is fermion number and the fact that the solitons have 
fermion number ±~ [|Tj| is the phenomenon of fractional fermion number. When A > 0, 
the ground state will fill with antisolitons. These scatter among themselves with the S- 
matrix element Q 

There is thus only one element in N a i, which is (t — 1) /2t. Using (|T7f), we see that 
e(— oo) = tA/(t + 1), and plugging this into ([TBI) gives the ground-state energy density at 
the critical point 

S(A - oo) = "^^ry = 

giving the correct result c = 1 and the correct j3 dependence of (|8|). This is thus a very non- 
trivial check on the sine-Gordon S-matrix and the relation (|H)p. In fact, this calculation 
can be alternatively viewed as a way of deriving the relation (|19|) between parameters t in 
the S'-matrix and /3 in the Lagrangian. 

Our next examples are the "fractional supersymmetric sine-Gordon" models @. The 
soliton content of these theories is described by k + 1 vacua on a line with "RSOS" solitons 
interpolating between adjacent vacua. For each pair of vacua ± 1) there is a doublet 
of solitons with charge ±i. The conjectured S'-matrix is a tensor product of the k — th 



RSOS [16] S-matrix governing scattering in the soliton labels and the ordinary sine-Gordon 



S-matrix governing scattering in the doublet labels H,|l5|. When A > 0, all the negatively 



charged particles can be in the sea. While scattering among these particles is diagonal 
in the ordinary sine-Gordon labels, it is not diagonal in the soliton labels. For example, 
the two-particle state with + 1) and (j + solitons can scatter into itself or into 
(j, j — 1) and (j — To quantize the momenta, one must find the analog of ( |10D in this 
non-diagonal case; in other words, one must find the eigenvalues for bringing one particle 
through an ensemble of the others. This is precisely what one does in a traditional Bethe 
ansatz calculation. In this case, the calculation has already been done in the minimal-model 



guise [p!7yiq|. As in the TBA, the effect of the soliton degrees of freedom is to introduce 
additional "pseudoparticles" . They enter into the analysis just like the real particles, but 
have zero mass and zero charge. The result is that there are k — 1 pseudoparticles (labelled 
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1 . . . k — 1) and one "real" particle with charge — ~ (labelled by k). This calculation gives 
the kernel <f> a ,a±i(@) = l/(27rcosh#) with all others zero; the only modification for our 
full tensor-product case is that there is an additional piece coming from the sine-Gordon 
S'-matrix. Since this part of the scattering is diagonal, (j>kk is given by (^0|). This yields 
N a ,a±i = \, — (t — l)/2£, and the rest are zero. The solution of ([17]) is 

atA 

'^-^ = wtty 

Only the real particle k has charge, so it is the only one which contributes to the energy 

(FBI) 



hi A 2 

E{A ->oo) = - . (21) 

When t=2 these theories are N=2 supersymmetric and thought to describe the 
integrable theories obtained by perturbing the k-th N=2 minimal model in its least relevant 
operator WB. Setting t=2 in ( PT| ) gives the correct central charge 3k /(k + 2). Comparing 
(H) with (|2l| ), we also find that the radius R of the bosonized U(l) current is related to 
the quantum-group parameter t by 

R 2 N=2 _ t(k + 2) 
R 2 ~2(t + k)' 

Setting R~ 2 =2 = 8ir(2/k + 2), this is in agreement with the relation derived using the 
quantum-group symmetry in (eqn. 3.12 there). The k—2 case corresponds to the iV=l 
sine-Gordon model. Another special case is the t — > oo limit, where the quantum-group 
symmetry becomes ordinary affine SU(2) k ; this corresponds to the SU(2)k WZW model 
perturbed by the current- current interaction. 

Another integrable perturbation of the N=2 discrete series is by the most relevant 
operator [j!9|,2~0|] . The vacua can be drawn as the k + 1 vertices of a regular polygon, labelled 
1 . . . k + 1 counterclockwise, and the spectrum consists of solitons interpolating between 
any pair of vacua. Each such soliton is a doublet; the charges of the doublet interpolating 
between the vacua j and (j + r) mod k + 1 are given by rj (k + 1) and r/(k + 1) — 1 |21[ 



With a positive magnetic field, one would expect that any negatively-charged particle (the 
member of each doublet with charge r/ (k + 1) — 1, which we can label by r = 1 . . . k) 
would be allowed in the sea. These scatter diagonally among themselves (so we don't need 
pseudoparticles) with the S'-matrix element a r ^ s {6) obtained in We need only the d rj i, 
which give 

(k + l-r)(k + 2) 



(k + iy 



from which it follows that the solution of fllTD is e r (— oo) = if r 7^ 1 and e\{— 00) 
(k + l)A/(k + 2). Using these in (0) gives 



^-») = -*T2S- (22) 

Only the r=l particle with the most-negative charge — /c/(/c + l) contributes to this leading 
term in the energy; in fact it can be shown that the other particles do not appear in the 
sea at all. 

Comparing (1213) with (^), we again obtain the correct central charge of 3k /{k + 2). 
This provides a check beyond the TBA [51] that the solitons form the entire spectrum of 
the theory, and the ^-matrix given there is the correct one. As for deforming the N=2 
symmetry into a quantum- group symmetry, this would presumably give an affine To da 
field theory with imaginary coupling constant and a background charge [ p0| , |2"2"|| , but these 
theories have not been studied a great deal. However, this calculation should provide the 
relation of the Toda coupling constant to the ^-matrix parameter. 

We also note that result (|j) can be continued to imaginary background field. Setting 
A = iaT and taking the temperature T — > 00, we reach the UV critical point. Multiplying 
the energy density from (|8|) by the box length 1/T (associated with finite temperature) to 
obtain an energy and adding the A=0 energy which follows from ||, we obtain the total 
ground state energy at the critical point 

This result is in perfect agreement with several specific examples [|ToT,|TB|| obtained by 



including imaginary chemical potentials in the TBA calculation and using dilogarithm 
identities. For a N=2 system, where (Rjy = 2/ R)=l, this is a familiar result for fermion 
boundary conditions twisted by e iaF . 

This of course does not exhaust the set of conjectured N=2 S'-matrices. We have 



checked that the conjectured S'-matrices for the N=2 sine-Gordon model |23[ also give 
the correct answer. There are models with more general soliton structure discussed in 



21|,E2| which have not been checked here or via the TBA. The background-field calculation 



described here will almost certainly turn out to be easier, because one only needs to 
diagonalize the S'-matrix for negatively-charged particles, not the full spectrum. In 
addition, there is no need for dilogarithm identities (although some may feel this is a 
disadvantage!). Another advantage of this method is that being at zero temperature, we 
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do not need to worry about finite-size effects. This is useful because there is at least one 
case where the presence of level crossings in finite volume means that calculating anything 
from the finite-size effects of the TBA is suspect ||. 
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